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ABSTRACT 



In the current work, the author present a symbolic algorithm for finding the determinant of any general 
nonsingular cyclic heptadiagonal matrices and inverse of anti-Cyclic heptadiagonal matrices. The algorithms are 
mainly based on the work presented in [A. A. KARAWIA, A New Algorithm for Inverting General Cyclic 
Heptadiagonal Matrices Recursively, arXiv:101 1.2306vl [cs.SC]]. The symbolic algorithms are suited for 
implementation using Computer Algebra Systems (CAS) such as MATLAB, MAPLE and MATHEMATICA. 
An illustrative example is given. 
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1. INTRODUCTION 

The n x n general periodic heptadiagonal matrices are takes the form: 
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where n > 8 . 

The inverses of cyclic heptadiagonal matrices are usually required in science and 
engineering applications, for more details (see special cases, [1-9]). The motivation of the 
current paper is to establish efficient algorithms for computing determinant and inverting 
cyclic heptadiagonal and anti-cyclic heptadiagonal matrices of the form (1.1) and (1.2). 

The paper is organized as follows. In Section 2, new symbolic computational algorithms, 
that will not break, is constructed. In Section 3, an illustrative example is given. Conclusions 
of the work are given in Section 4. 

2. Main results 

In this section we shall focus on the construction of new symbolic computational 
algorithms for computing the determinant and the inverse of general cyclic heptadiagonal 
matrices. The solution of cyclic heptadiagonal linear systems of the form (1.2) will be taken 
into account. 

Throughout this section, the parameter t is just a symbolic name and detH is the 
determinant of the heptadiagonal matrix of the form (1.1). We state the following result 
without proof (see [10,1 1]). 

Theorem 2.1. Suppose tha: 
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Then, we have the two-term recurrence 

u i =a i u i _ 1 , i =1,2, ...,n 
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At this point it is convenient to formulate our first result. It is a symbolic algorithm for 
computing the determinant of a cyclic heptadiagonal matrix H of the form (1.1) 

Algorithm 2.1. To compute det H for the cyclic heptadiagonal matrix H in (1.1), we may proceed as 
follows: 

Step 1: Set ai=di. If c<i=0 then cti= t end if. Set Ui=di, gi=ai, Zi=Ai, k^An.i/ di, Vi=bi, Wi=Bl, h!=a n / 
ai, Wi=B! , f 2 =b 2 / ai, e 3 =B 3 / ai, a 2 =d 2 -f2 *gi- If a 2 =0 then a 2 = t end if. Set K 2 =-ki*g!/ a 2 , u 2 = 
a 2 *u 1 ,v 2 =B 2 -f 2 *v 1 , w 2 =-f 2 *Wi, h 2 =(A n -hi*gi)/ a 2 , g 2 =a 2 -f 2 *Z!, f 3 =(b 3 -e 3 *gi)/a 2 , a 3 =d 3 -e 3 *z 1 -f 3 *g 2 . 
If a 3 =0 then a 3 = t end if. Set u 3 = a 3 *u 2 , k3=-(k 1 *z 1 +k 2 *g 2 )/a 3 , h^^hxZi+h^)/^, v 3 =-e 3 *v r 
f 3 *v 2 , \N i =-f i \N 2 -e i \N 1 . 
Step 2: Compute and simplify: 
For i from 4 to n-2 do 
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Step 3: Compute and simplify: 
For i from 4 to n-5 do 

K i =-(k, 3 *C, 3 +k i . 2 *z, 2 +k i . 1 *g, 1 )/ a, 
Wi=-(Di*Wi_ 3 / a^+e^w^+f^WM) 
End do 

Step 4: Compute and simplify: 
For i from 4 to n-4 do 

h i =-(h i . 3 *C i . 3 +h i . 2 *z i . 2 +h i . 1 *g i _ 1 )/ a, 
Vi=-(Di*Vi_ 3 / a^+e^v^+f^v^) 
End do 

Step 5: Compute simplify: 
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hn- 3 = (D n -h n _4*g n _ 4 -h n _5*z n _5-h n _6*C n _ 6 )/ a n - 3 
h n - 2 =( B n -h n . 3 *g n . 3 -h n _ 4 * z n . 4 -h n _ 5 * C n _ 5 )/ a n _ 2 
v n _ 3 =C n _ 3 -D n _ 3 *v n _ 6 /a n _ 6 -e n _ 3 v n _5-f n _ 3 v n _ 4 
Vn- 2 =A n _ 2 -D n _ 2 *v n _ 5 /a n _5-e n _ 2 v n _4-f n _ 2 v n _ 3 
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Step 6: Compute det H =(u n ) t = . 

Now we formulate a second result. Following [10], we can using CHINV algorithm to 
compute the inverse of a general cyclic heptadiagonal matrix of the form (1.1) when it 
exists. The following theorem given the inverse of cyclic anti-heptadiagonal matrix H antl . 



Theorem 2.2. Let P be an n X n matrix as the following form 
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Then the inverse of cyclic anti-heptadiagonal matrix H antl is given by 

j_£ and _ pj_£ 1 

3. An illustrative example 

In this section we give an example for the sake of illustration. 
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Example 3.1. Consider the 10 X 10 cyclic heptadiagonal and cyclic anti- heptadiagonal 
matrices 
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By using algorithm 2.1 and theorem 2.1 we obtained the determinant of matrix (3.1) and 
inverse of matrix (3.2) respectively: 
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4. Conclusions 



In this work new computational algorithms have been developed for computing the 
determinant and inverse of general cyclic anti-heptadiagonal matrices. The algorithms are reliable, 



computationally efficient and will not fail. The algorithms are natural generalizations of some 
algorithms in current use. 
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